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The Landauer's formula conductance of the telescoped 
armchair nanotubes is calculated with the Hamiltonian de- 
fined by first-principles calculations (SIESTA code). Herein, 
partially extracting the inner tube from the outer tube is 
called 'telescoping'. It shows a rapid oscillation superposed 
on a slow oscillation as a function of discrete overlap length 
fl/ — i j a with an integer variable L and the lattice constant 
a. Considering the inter layer Hamiltonian as a perturbation, 
we obtain the approximate formula of the amplitude of the 
slow oscillation as |j4|^/(|ylp -\- e^) where A is the effective 
interlayer interaction and e is the band split without inter- 
layer interaction. The approximate formula is related to the 
Thouless number of the dispersion lines. 



I. INTRODUCTION 

Single-wall nanotubes (SWNT) show various useful 
characteristics as nanoelectronic devices. ^'^ They can be 
either metallic or semiconducting, depending on their 
chirality.^ Although a semiconducting SWNT can act as 
a transistor by itself,'* we can increase NT device func- 
tionality by assembling SWNTs. For example, a junc- 
tion between a metallic SWNT and a semiconducting 
SWNT can work as a Schottky diode.^^^ The SWNTs 
can be connected by covalent bonds^'^ and interlayer 
interactions. ^^^ Double- wall nanotubes (DWNTs), multi- 
wall nanotubes (MWNTs), and nanotube bundles are 
all examples of SWNTs that can be naturally assem- 
bled through interlayer interactions. Processing these 
naturally assembled SWNTs is one way of constructing 
higher-order structures for NT devices. One prototype 
higher-order structure is the telescoped double wall nan- 
otube (TDWNT), in which the inner tube is partially 
extracted from a DWNT, as shown in Fig. 1. Since the 
interlayer force is relatively small, the inner tube can be 
slid and rotated with respect to the outer tube by at- 
tached piezoelectric probes with a negligible change in 
intralayer bonding. ^'^^ In this way, we can realize vari- 
ous interlayer configurations in the TDWNT, unlike the 
case for bulk graphite. When a bias voltage is applied 
between the attached probes, the resulting current is sig- 
nificantly influenced by the interlayer configuration be- 
cause a carrier from one probe must travel an interlayer 
path to reach the other probe. This sensitivity can be 
used in nanomechanical switches and nanodisplacement 
sensors. In addition, TDWNTs enable us to measure the 
interlayer conductance more directly than un-telescoped 
DWNTs." 

Among the reported theo- 



retical work on TDWNTs, *^~*^, there are controversial 
results regarding the conductance of a TDWNT com- 
posed of (10,10) and (5,5) armchair nanotubes. Com- 
pared to a quantum conductance Go = 2e^/ft., the max- 
imum conductance was found to be 2Gq based on tight 
binding (TB) model calculations, but only Go based on 
first- principles calculations.*'*'*'^ Landauer's formula was 
used in these calculations, so the conductance in units of 
Go equals the sum of the transmission rates. In Ref.*^, 
the author derived approximate formulas to clarify the 
direct relationship between the transmission rate and the 
interlayer Hamiltonian. The relation between the trans- 
mission rates and the dispersion lines was used in this 
derivation. The energy E, Hamiltonian H , and overlap 
matrix B were common to the exact calculations and the 
approximate formulas. The approximate formulas indi- 
cated that the disagreement comes from differences in 
the interlayer Hamiltonian. The interlayer Hamiltonian 
is also important in multilayer graphenes*^ with regard 
to, e.g., transmission through a boundary between mono- 
layer and bilayer graphene.*^ The TB Hamiltonian was 
used in Ref.*^, but is insufficient for the interpretation 
of the first-principles calculations. Thus, both the exact 
method and the approximate formulas were adapted to 
a first-principles calculation code with an atomic orbital 
(AO) basis. The effectiveness of the approximate formu- 
las was judged from their agreement with exact calcula- 
tions. 



II. INITIAL GEOMETRIC STRUCTURE AND 
SITE INDEX 

Geometric optimization was performed with the ini- 
tial structure defined as follows. In the initial structure, 
cyhndrical coordinates {r,9,z) — {rf^ ,0'im^^'j,m) of the 
carbon atoms in the tube fj, — O (outer tube) and fJ- = I 
(inner tube) are represented by 
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The notation of Eqs. (1) and (2) is the same as Ref.*^ 
except that the integer index j in Ref.*^ is replaced by 
2j -|- m here. The radius of the tube fi is denoted as 
r^ — v3an^/(27r), with integers no = 10, n/ = 5 and the 
lattice constant a = \/3 x 0.142 nni. The unit cells are 



numbered j along the tube axis, while the atoms in each 
unit cell are numbered m = —1,0 and I = 1, 2, • • • , 2n^. 
As the DWNT considered here has no mirror plane par- 
allel to the tube axis, symmetric and anti-symmetric 
states are hybridized by the interlayer interaction. We 
use /i, l,j,m as site indexes in the following. We discuss 
three structures: an un-telescoped DWNT with five unit 
cells , a TDWNT, and an un-telescoped DWNT with in- 
finite length. These are denoted as [fi], [td], and [in]. 
Although Eqs. (1) and (2) are common among [fi], [td], 
and [in] , the range of j at which the corresponding atoms 
exist is different, as 



1 < J < 5 



[fi]. 



and 



1 < j < oo in tube O of [td] 
-oo < j < L in tube / of [td] 



-oo < J < oo m [mj 



(3) 



(4) 



(5) 



where L denotes the number of unit cells in the DWNT 
of [td] . In the initial structure of [fi] , the direction of the 
C-H bond is the same as the corresponding C-C bond 
and the C-H bond length is 0.11 nm. 



Eq. (6). The ranges of j' and j in Eq. (7) are shown 
in Eqs. (4) and (5). From the obtained H^'ft '' ' and 
B'^'!^ , the block matrices of structure [in] are defined 



as 






tTjt^ ,M 
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H: 



3,3+Aj 



ioI{^l,^i') = {IJ),{o,I),{I,o), 



h: 
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H. 



0,OJfi] 
3,3+Aj 



AsB. 



0,OJfi] 
3,3+Aj 



and 






^ri+Aj,rj 



^3,3+Aj 



(8) 



(9) 



(10) 



where i] is an arbitrary integer and Aj = 0, 1,2. When 



|Aj| > 2, Bi;^^^^ 



and H: 



Ai'.AiJ?] 



0. Here we select 



'7,'7-l-Aj 

block matrices H^^\B^^^ for which the boundary effect 
is relatively small and build them into the H and B ma- 
trixes of structure [in]. Once we obtain H^^^ and B^^^ 
from the SIESTA, there is no adjustable parameter ex- 
cept for Ae. The purpose of introducing the parameter 
Ae will be explained in Sec. IV. The H and B matrices 
of structure [td] are defined also by the building block 
procedure as will be shown in Sec.VB. 



III. BUILDING BLOCK PROCEDURE 



IV. DISPERSION RELATION 



The SIESTA code is used as the first-principles calcu- 
lation of structure [fi].^*"^" The geometry is optimized 
by the conjugate gradients method with the convergence 
criteria 0.04 eV/Ang. The optimized and initial struc- 
tures are superimposed on Fig. 2. We can see that the 
initial structure is almost maintained after the optimiza- 
tion. The minimum, average and maximum bond lengths 
are 0.1175, 0.1177, 0.1181 nm for C-H and 0.1413, 0.1451, 
0.1482 for C-C. The output data after the geometrical 
optimization are used. Since the obtained Hamiltonian 
matrix has an atomic orbital (AO) basis, it can be di- 
vided into block matrices based on the unit cells. The 
one-electron wave function -0 is represented by the linear 
combination of AGs {0}, as 

^= E EE E E<™.«.<^"™,.. • (6) 

fj,—I,0 j l—l m— — 1 a 

In Eq. (6), four single zeta AOs {a — 2s, 2px, 2py, 2pz) 
are considered per atom. 

The secular equations for structures [in] and [td] are 
represented by 

E E^/^'-^ E E^^/'^' (7) 

fi'=I,0 j' fi'=I.O j' 

where H'^'P and B^'v are the block matrices of H and 
B, respectively. The vector d -^ is composed of dj'^^^^ ^ in 



In order to discuss the dispersion relation of structure 
[in], it is useful to define 



X 



D.D _ 



and 



XOJ 



M 



XO.o 



X^'-^ik)^ J2 <, 



fJ. ikaj 



(11) 



(12) 



i=-M 



where X = H,B. The dispersion relation El^{k) and the 
wave function per unit cell g!^{k) can be obtained by 



H^^^{k)rAk) - EAk)B^'Ak)gAk) 



(13) 



with the band index r. According to Eqs. (11), (12) and 
(13), E!^{k) corresponds to a DWNT with interlayer in- 
teraction, while E^k) and E^{k) correspond to isolated 
SWNTs. 

The dispersion lines near the neutral Fermi level Ep 
are denoted by E^^, E^^, E^3, £^±4 for a DWNT with 
interlayer interaction, and E^i, E^i E^2 ^'^d ^±2 lor 
isolated tubes O and /. Figure 3 shows the dispersion 
lines numbered based on the index r under the conditions 

K\ ((-i)^t) - ^-M (-(-i)^if ) - 4°^ -d r^ > 

0. For the isolated tubes fi — I and fi — O, Ef^^ and E^2 
correspond to the symmetric and anti-symmetric states. 



respectively, with respect to the mirror plane parallel to 
the tube axis. For a DWNT with interlayer interaction, 
however, the mirror symmetry is broken and the band 
splitting is increased as E^ - E[' > \E^ - E(\ and E^ - 
E^>\EO-Ei\. 

Figure 4 shows the dispersion relation of the isolated 
SWNTs for M = 1 (open circles) and M = 2 (closed 
diamonds). The shift between M = 1 and M = 2 is 
about 0.002 - 0.003 Ry for the E^i band. In order 
to include this energy difference, the integer M in Eq. 
(12) is set to 2 hereafter. Squares and circles in Fig. 5 
show the intralayer and interlayer effective Hamiltonian 
elements h{r,r') , respectively, for the atomic distance 
0.3 nni < |r — r'\ < 0.45nm. The effective Hamiltonian 
h{f, r') is defined by 



l^a.a' 3r,a-"5r',Q ' 



l^a.a' 9f^a^9r,a' J Z^a,a' 9p'^a^9r',a' 



(14) 



where gp,a denotes the component of gLi,g^i for 
k — 0.667r/a, atomic position r and AO type 
a=2s,2pa;,2pa;,2p2. The spatial range of the effective 
Hamiltonian is considerably longer compared to the TB 
Hamiltonian used in Ref.^^. Crosses represent the in- 
tralayer elements between the second nearest neighbor 
unit cells. They cause the energy difference between 
M = 1 and Af = 2 in Fig. 4. 

According to Ref.^'"', we define the intrinsic shift as 



er^E?{k)-E',{k) 



(15) 



Figure 4 shows a reasonable value of ^r- However, 
|er(fc)| is considerably larger than other first-principles 
calculation results, probably because of the finite size ef- 
fect of structure [fi] used in the first step of the building 
block procedure. The parameter Ae in Eq. (8) changes 



as 



^,(0) 



Ae 



(16) 



where e)- denotes the intrinsic shift in Fig. 4. As Ae is 
common to the exact and approximate methods, it is not 
a parameter that can be adjusted to produce agreement 
between the results obtained by the two methods. Figure 
4 shows that eL°l = -0.022 Ry and £^°^ = -0.028 Ry at 
k = 0.657r/a. 

When ^ = D, Eq. (13) can be transformed to 

V.^(fc)* [i?^'^(fc) - E?,{k)B^^^{k)] g^,{k) = . (17) 

In the following, we consider the wavenumber fc as a func- 
tion of the energy E. The wavenumbers of the isolate 
tubes / and O are denoted by k^ {E) and k'-^{E), respec- 
tively. Neglecting mixing between different r values, the 
lowest order approximation is represented by 



'gJ^ 



{k)c,{yj'g^{k'),yo'g^''{kO)) 



(18) 



where r' — r, r + 2 -pr . Here the three different wavenum- 
bers have the common energy as _E = E^(k) = E^{k^) — 
E'='{k'=') and k^y ^ fcO. 
Using Eq. (18) and 

X^''■^'{k)~X^''■^'{k^') + {k-k^')^X^''^^' , (19) 

Eq. (17) can be approximated by 

Rc[{A^ + C.)y*oyi]+ E ^5^^|y,p^0 



where 



5^ = *5^(fc^)*B^^^(fc^).g('(fc^) . 



b't 



dEl _ ( t^^.. dP^^^{k) ^^ 



9' 



9^ 



k=kf 



Ar 



dk \ "^ dk 

2'g^O{kO)*pOJ(k^)g^i(k^) 



VWK 



a = 



- 2(^ ^') *5.°*(fc?) ' "^ -^ 



dpo^ 



dk 



9r 



VWK 

P^'^'(fc) = Hf^'f^'ik) - EBt^'f^'ik) 



(fpM,M (iiyp.p dB^'^^' 

— E- 



k=k' 



dk 



dk 



dk 



Here Eq. (22) will be proved in Appendix B. 
From Eq. (20), k^, can be approximated by 

k^^^^iE) ^ \ [ki{E) + k?{E) ± Afcf '^] 



where 



Ak^^^ = \Akr\VTT^ 



Akr = k^{E) - ki{E) 
and Xt- is a dimensionlcss parameter defined by 
dk^ dk'=' 



Xf ^ -^7 



dE dE 



{Aky 



(20) 
(21) 
(22) 

(23) 

(24) 

(25) 
(26) 

(27) 

(28) 
(29) 

(30) 



The band split along the k axis without interlayer inter- 
action is denoted by Eq.(29), while that with interlayer 
interaction 



Afcf 



k^^,.iE)-ki^{E) 



(31) 



is approximated by Eq.(28). Since \C\ <C \A\, the effect 
of Eq.(24) is neglected here. Note that we define g^ , g^, 
V^ h\. , k^ and k\ without the interlayer matrixes H ' 
and B^-f. 



V. TRANSMISSION RATE 

In the following, T^-^t' denotes the transmission rate 
from the t'(= 1,2) channel of tube / to the t{— 1,2) 
channel of tube O. Here, the notation of r is common 
with that of Sec. IV, i.e., t = 1 and t = 2 are symmetric 
and anti-symmetric channels, respectively. 

Since the interlayer Hamiltonian breaks the mirror 
symmetry, the inter-channel transmission rates Ti,2 and 
'?2,i are not zero. Nevertheless, the large difference be- 
tween fci ~ — 27r/(3a) and ^2 — 27r/(3a) suppresses T12 
and T2.1. In the approximate formulas, we consider only 
the intra-channel transmission rates Ti 1 and T2.2. 



A. Approximate formula 



is proportional to the geometrical mean of the density of 
states. Since the dispersion lines are almost straight and 
parallel. 



dk 
dE 



Akr 



(37) 



Using Eqs. (28), (30), (36) and (37), we can transform 
Eq.(35) to 



T' 

T.T 



Ak-r 



Ak?'^ 



sm La 



V 2 



1 + x, 



■sin \/l + Xr — ^La 



(38) 



In order to clarify the relation of Eq. (23) to the ef- 
fective interlayer interaction of Ref.^^, we transform Eq. 
(23) into 



A = 



2jx^{lT{H-E)'^'{k)d^f 



^J xO(l)**0(fc)dVy^/ x-f (l)**^(fc)dV 

where H is the Hamiltonian operator, 

00 



(32) 



j = -oo 



2n„ 



x"(j)^E E E<™.«,.w<™,a,. 



(33) 



(34) 



/=1 m=— 1 a 



and g^^ ^ Ak) denotes the component of g^{k). Equa- 
tion (32) clearly indicates that ^ equals the matrix ele- 
ment of the perturbation H"^'^ — EB'^'^ between the un- 
perturbed states ^*^ and '^^ per unit cell. When jj, — 0,I, 
Eq. (33) is the eigenfunction of the intralayer Hamilto- 
nian iJ^''^ and bears no relation to the interlayer Hamil- 
tonian H'^-^ . Since A of Ref.^^ can also be represented 
by Eq. (32), the physical meaning of A is common to 
the present paper and Ref.^^. There is, however, another 
possible definition of A where P'-''^ = H'-''^ — EB^-^ is 
replaced by H '^ . The two definitions cause no differ- 
ence when the interlayer overlap matrix B"^'^ is absent. 
Though the effect of EB'~''^ is also investigated in Sec. VI, 
the definition (23) is chosen here. 

The approximate formula of Ref.^^ is represented by 



T', 



where 



l^rl 



\A, 



dk 

'dE 



■dk 



V/^^^W^^2 



dk? dki 



dE dE 



(35) 



(36) 



B. Exact numerical calculations 

The scattering matrix method was used as an ex- 
act method. ^^~^^ We can also use the Green's function 
method, which is useful for including inelastic scattering, 
electron correlation and finite bias.^^'^^ Without these 
effects, however, the same transmission rate is obtained 
by both methods. ^^ There are advantages of using scat- 
tering matrix method compared to the Green's function 
method. One is the explicit relation to the wave func- 
tion, and the other is that we can estimate the numerical 
error from the unitarity of the scattering matrix. 

We can obtain only the propagating states from Eq. 
(13), but the evanescent wave states are also necessary for 
the exact calculation of the transmission rate. In order to 
obtain both the propagating and evanescent states, the 
transfer matrix T ^ is defined as 



where 



Tm^ 



n 



3 



/ 0, 1, 0, \ 

0, 0, 1, 

0, 0, 0, 1 

I Vt^ v^i- v^i- vt^ , 

\ ^2 ' ^i ^ -^0 ' -^1 / 



(^r) 



im) 



pp. A* 
3,3+i 



■^3,3+j 



The matrix P is defined as 



pM'.P — rr^'./J Jfi] 



1,1 



f ,1 



{E + Ae5 



fj,%iiSfj,,o)B':^, 



/j\/jjfi] 



(39) 



(40) 



(41) 



(42) 



Figure 6 shows the relation between structure [fi] and 
the matrixes P. The matrix P^-^ is composed of 

plj^ pOJ^ phO ^-^^ pO^O ^^ ^Jjg gj^j^g ^g^y g^g ^^ ^^l). 

The secular equation is represented by e^+i = F^ej 
where *ej = {^df^^i *'^i'^ii ^d.^^, ^df^_^^) with /i correspond- 
ing to j as (^ = I,j < -2), in ^ D,3 < j < L-2) 



and {jj, = 0,L + 3 < j). For the propagating waves, 
the eigenvector u and the eigenvalue A of the trans- 
fer matrix F^ satisfying F^u = Xil and Cj — \^u can 
be related to k and 5^ of Eq. (13) as A = e**^" and 
*u = (t^A'^ A*5^, A2*(7^, A^V^). In the following, the 
dimension of g^ is denoted by iV^ {Ni ~ 80, Nq — 
160, Nd = 240) . The 4iV^ independent eigenvectors of 
4iVp X 4iV^ matrix F^ are labeled by r as follows. For 
the propagating waves {\X'^\ = 1 ), t = ±1,±2, • • • ,±NI^. 
For the evanescent waves (|A'5^| 7^ 1 ), r = ±(^^ + 
l),±(iV^ + 2),---,±27V^. The sign of t is chosen to be 
consistent with the propagation direction and decay di- 
rection when \t\ < N^^ and \t\ > N^^ + 1, respectively. 
Here, the channel number A^^ denotes half the number 
of independent propagating waves in region /i. In the 
present paper, the energy E is close to the neutral Fermi 
level, and thus N^ = iVf = 2 and A^^^ == 4. 
The amplitude vector dj is represented by 



Wu 






J2 (Vr5r"7r" 



t^\^ nt^rv M 



(43) 



-—Na 



where r 7^ and the correspondence between fi and j is 
represented by (p = I,j < 0), {n = D,l < j < L) and 
(/i = 0,L+l<j). 

In the periodic region (,u = I,j < —2), {^ ~ D,3 < 
j < L — 2) and {n — 0,L + 3 < j), the secular equations 
can be represented by 



^M - \^ Y^ d 






E ^A^.^; 



J+Aj 



Aj=0 



with Eqs. (40) and (41). In the transition region, the 
secular equations can be represented by 



*Q(ci3 + O^cii + qiD^i'' = 



and the matrixes q are defined as 



% 



pD.D 

tpD,D 
-'^1.2 ' 



pD,D 
-^1.2 
pD,D 
-'^2.2 



g/D 




II 

% 



pD.D „D,D 
-^4.5 ' -^5.5 



tr>I,0 tr>0,0 tpI-.O tr>0,0 



*P. 



QDO 



3,5 





'P' 



3,5 





p 

3,5 ' 



*P. 



3,5 



(52) 



(53) 



(54) 



(55) 



Although the coefficients {7^}, {7^} and {7^} in Eq. 
(43) can be chosen arbitrarily to satisfy Eq. (44), they are 
subject to the conditions expressed in Eqs. (45), (46), (47) 
and ,(48). Using Eqs. (43), (45), (46) and (49), we can 
obtain a 320 x 320 matrix a^^ satisfying CjI,*^+) = 
(*7/,*7-^)*'T^^, where*7^ ^ (7^1,7^2, ••• ,7^^,)- The 
matrix ct^° satisfying (*7_^,*7_f) = (*7^,* 7°) *cr^° 
can be obtained in the same way. By eliminating 7_|? 
from these equations, we can obtain the matrix a^ satis- 
fying (*7i,* 7+') = (*7-/,* j_^Ya'^. The 4x4 scattering 
matrix S^^ is obtained from the 240 x 240 matrix cr^*^ 



(44) as {S'O^r.r' = (^^°)r,r', 8^%,^, = a 



- ^lO 



'r+80,T' + 80 



lO qlO _ 

r+80,T'' '^T,r'+2 ^ 

, where r = 1 , 2 and 
t' = 1,2. From the scattering matrix, we can obtain 

Tt.t' = 15^+2,1-'! • 

The difference between ^S^'^*S^'^ and the unit matrix 
(45) is represented by 



''t,t'+80 ^^'i '^t+2,t' 



*9/Z5cii + g^ci^ + Qf C3^ = 



(46) err(5^0) ^ ^ E 



1^=112 = 1 



V^ qIO* qIO 



'^hM 



(56) 



'Q?cE-z + %cE-i + qoocE+i = (47) 



^qoocE-i 









where the matrixes Q and the vector c are defined as 



^3 — 



^3^ 



(49) 



When there is no numerical error, 5*^*^ must be uni- 
tary and Eq. (56) must be zero with the normalization 
shown in Appendix A.^^'^^ As the unitarity represents 
that the sum of transmission rate and reflection rate must 
be unity, the numerical error in the transmission rate is 
estimated to be Eq. (56). 



VI. RESULTS 



r,^ 



Q', 



pM,M 

-^3,3 ' 

t pfJ-.t^ 

-'^3.4 ' 



pt^,t^ 
-^3.4 
ptJ-^t^ 
-^3.3 



pM,M 

pM,M 
-f^3,4 ' 



pM,M 
^3,5 



Figure 7 shows the dispersion lines for (a) Ae = —0.015 

(50) , (b) Ae = 0, and (c) Ae = 0.020 Ry. The relationship be- 
tween Ae and Sr is shown by Eq. (16). At fc = 0.657r/a, 
e_i = -0.022 + Ae and e2 = -0.028 + Ae. Crosses repre- 
sent E^{k) and E^{k) without interlayer interaction. In 

(51) pjg Y ^^^^ jt;0_^ and ELi almost coincide with each other 



as £_i ~ 0. The dispersion lines with interlayer inter- 
action are shown by circles for the exact E (k) and by 
solid lines for the approximate E^'^{k) defined by Eqs. 
(27), (28), (29) and (30). Since jAal < |£2|, E^ and E^ 
are almost the same as E2 ' and E^ ' , respectively. In 
the symmetric bands specified by t — —1, on the other 
hand, the relatively large |A_i| causes a clear difl^erence 
between the solid lines and the crosses. As |A_i| is al- 
most constant, E'^-^ and E^^ are almost parallel to E_{ 

and E_^ . The solid lines reproduce the band split of 
T — —1 channel, but underestimate that of r = 2 chan- 
nel. 

Figures 8 and 9 shows Tii and T2.2, respectively, as a 
function of the number of unit cells L. The circles rep- 
resent the exact transmission rate. The values of Ae in 
Figs. 8 and 9 are the same as in Fig. 7. Because L 
must be larger than three in the exact calculation with 
Eqs. (46) and (47), i > 4 in Figs. 8 and 9. Triangles 
in Fig. 9 represent Eq. (56) showing the numerical er- 
rors of Til + 72.2 + ^2^1 + ri.2 in the exact calculation. 
The average and maximum of Eq. (56) arc 4.9 x 10""* 
and 1.4 x 10~^. The magnitude of T1.2 + r2,i obtained 
from the numerical calculation is also 10^^ ^ 10^*^. This 
indicates that Ti_2 + ^2^1 cannot be discussed with suffi- 
cient accuracy. Thus we focus our discussion on Ti^i and 
T2,2. Although the interlayer Hamiltonian could cause a 
pseudo-gap in the energy band of the DWNT,^* there is 
no pseudo-gap at the energy E — —0.4 Ry chosen here. 
The circles shows a rapid oscillation superposed on a slow 
oscillation. The rapid oscillation is due to a standing 
wave caused by a large intralayer reflection at the open 
edges j = l,i in structure [td]. As the wavenumber k is 
close to ±27r/(3a), the period of the rapid osciUation is 
close to three. In order to see the slow oscillation of the 
exact r(L), 

T"-(L) ^ nL-l)+T{L)+TiL + l) ^^^^ 

is shown by the dot-dashed lines. In order to estimate the 
edge effect, q'^ in Eq. (46) and Qq in Eq. (47) are replaced 
by Qq . The exact Ti^i and T^Y with this replacement 
are shown by crosses and dashed lines, respectively, in 
Fig. 8. Agreement between the dashed line and dot- 
dashed line indicates that the slow oscillation defined by 
(57) is not sensitive to the detail of the edge modeling. 
Though crosses and dashed lines are omitted in Fig. 9, it 
is also confirmed that the two edge models show similar 

rpave 
^2,2 • 

The thick solid lines represent T' defined by Eq. (38), 
while the thin solid lines represent Eq.(38) where EB'-''^ 
in Eq.(23) is replaced by zero, i.e., the interlayer overlap 
matrix is set to zero. The thin solid lines are shown only 
for the first period. Though the thin solid lines are close 
to T"''"^ for the small L, the thick solid lines qualitatively 
reproduce T'""^ for the wider range of L. Thus we con- 
centrate our discussion to the thick solid lines below. The 
effective interlayer interaction of the symmetric state Ai 



is much larger than that of the anti-symmetric state A2. 
This is the reason why Tii^T2.2- In the derivation of 
Eq. (35) shown in Ref.^^, the interlayer Hamiltonian is 
the same as that for structure [td], while the intralayer 
Hamiltonian is the same as that for structure [in]. Thus, 
Eq. (35) does not include the strong intralayer reflection 
at the edge j = 1,L of structure [td]. This is why the 
rapid oscillation does not appear in Eq. (38). 

The period of the exact T""^ is close to 27^/ /S.k^ where 
Afe^ denotes the exact band split of channel r with in- 
terlayer interaction. Since Afci ' ~ A/cf , the thick solid 
line and the dot-dashed line have almost the same period 
in Fig. 8. On the other hand, the thick solid lines show 
the longer period than the dot-dashed lines in Fig. 9 be- 
cause Afcg ' < Afc^. Nevertheless the peak height of 
rpave jg qualitatively reproduced by x/{l + x) of Eq. (38) 
both in Fig. 8 and in Fig. 9. 



VII. DISCUSSION 

In the present paper, we consider the differential con- 
ductance dl/dV at zero bias voltage V = 0. Wc have 
to use more sophisticated numerical methods to obtain 
the effect of the finite V on the current I.^^ If the rigid 
band picture is effective, however, the current / for finite 
V can be approximated by 



/ = 



h 



Ep+eV 



T{E,e + eV)dE 



(58) 



Ef 



where Ep is the Fermi level of tube O and the interlayer 
transmission rate is represented by T(E, s) as a function 
of the energy E and the intrinsic shift e. 

The effective interlayer interaction A defined by Eqs. 
(23) enables us to analyze the effects of the interlayer 
Hamiltonian on the transmission rate. For example, we 
can discuss the difference between the TB used in Ref.^^ 
and the SIESTA used here. The longer cut-off distance 
of the interlayer Hamiltonian decreases |yl2| as was dis- 
cussed in Ref.^^. The cut-off distance for SIESTA is 
much larger than that for TB as seen in Fig. 5. This 
is why IA2/A1I and r2,2 for SIESTA are negligible com- 
pared to those for TB. Using the perturbation formulas, 
we can distinguish the effects of the intralayer Hamilto- 
nian H^'^ , H^'^ from those of the interlayer Hamiltonian 
H'^'^ . Rigorously speaking, the interlayer interaction in- 
fluences H^'^ and iJ*^'*^ in the self consistent calculation 
of structure [fl] for the building block procedure. Never- 
theless, H^'^ and H'^'^ are approximately considered as 
the Hamiltonian of isolated SWNTs. As long as a real 
space basis such as an AO basis is used, the calculation 
methods in the present paper can be applied to various 
other first-principles calculation codes. 

Equation (38) is effective both in the SIESTA Hamil- 
tonian and in the TB Hamiltonian. We have to be care- 
ful that the effective interlayer interaction (23) must in- 
clude the interlayer overlap matrix B^'^ that is absent in 



Ref.^^. Using Eq.(38), we can predict the amplitude of 
T"''"'^{L) as seen in Figs. 8 and 9. Though the precision of 
Eq.(38) is lower in Fig. 9, the qualitative dependence on 
e is reproduced. The period of Tf^^(L) is also reproduced 
well, but that of T2Y is overestimated. It is common to 
TB and SIESTA tliat Eq.(38) is more effective for Ti^i 
than for T2,2- The period 27r/Afc has been discussed in 
other papers, but the amplitude formula a;T/(l + a;,-) in 
Eq.(38) is proposed firstly here. Note that the dimen- 
sionless parameter Xr is defined by Eq. (30) without a 
fitting parameter. Owing to the linear dispersion, the 
amplitude of T[ ^ are represented by the energy ratios as 



Afu Afu 



Aki 



Afcf 



£1 



Sf - E^ 



(59) 



When |£i|/(-Bf 



EH 



HEP 



close to 
Ell is the 



: 1 , Eq.(59) 
.3 -i^f-|ei|)/|ei|. Since E^-E[> 
energy shift caused by the interlayer interaction, its ra- 
tio to |£i| can be considered as the Thouless number. ^^ 
Author expects that Eq. (38) represent the generalized 
Thouless number analysis and can be applied to other 
related systems with commensurate interlayer configura- 
tions. 



APPENDIX A: RELATION BETWEEN THE 
NORMALIZATION AND PROBABILITY FLOW 

The probability flow J(7'|7) from orbital 7 = 
(/i, /,m, a,j) to orbital 7' — {jjl' ,1' ,ra' ^a' ,j') is repre- 
sented by 



J(7'|7) =; Im \{H^'^~f - EB^i^y)d*^,d^'\ 



(Al) 



where H, B and E represent the Hamiltonian matrix, the 
overlap matrix and the energy, respectively. Equation 
(Al) satisfies the two necessary conditions for station- 
ary flow: 'direction of flow' represented by 0/^(7] 7') = 
~J{l'\l) ^'^d 'conservation of probability ' represented 

by 



E^(Vl7) = 



(A2) 



Equation (A2) is derived from the secular Equation 
J^ji^-y'.-y ^ £^-67', 7)^7 = 0. Using Eq. (Al), the proba- 
bility flow Itotij) through the cross section a.t z — {j — ^)a 
is represented by 



ItotU) =Ij,j-l 






(A3) 



where 



^L' - Z! Z! J{tx,l,m,a,j\fiJ',m',a',j'). (A4) 

Corresponding to M = 2 in Eq. (12), 7^-/ = when 
\j~j'\ > 3. Using Eqs. (42), (43), (51), (Al) and (A4), 
Eq. (A3) can be expanded as 



iMj)^ E E 1- 



Sr.r' 



a , Li 



T = — Nfi t'=t 

where t 7^ 0, r' 7^ 0, 

r;,^^,^Im[(A^rAi!)-''z,^,,,7.";7(' 



K'.r ^ wr 



-2 (t r'fj.* 



*Q^// 



and 



//^ 



K9r 



(A5) 



(A6) 



(A7) 



(A8) 



In Eq. (A5), the correspondence between fj, and j is 
represented by (/z = I,j < —2), {fJ, = D,3 < j < L — 2) 
and (p = 0,L + 3 < j). 

On the other hand, we can obtain 

((Ai!)-2*Q5' + O(; + (Ai!)2Q5')./t = (A9) 

from the secular equations. Multiplying Eq. (A9) by 
*/t'^*, we can obtain 

(A^A^r)' z!;*^, + z^,^^ + *./vr o^/> = o . (aio) 

Exchanging t and t' in the complex conjugate of Eq. 
(AIO), we can also obtain 

(A^A^r)' z!;,^^ + z^;^, + *./vr o^/> = o . (ah) 

When the condition 

(A^A^D* = 1 (A12) 

does not hold, wc can show 



T ,T T.T 



[rrKrf + i 



(A13) 



from Eqs. (AlO) and (All) . Thus, ■q'^, .^ + ri^.^, has a 
nonzero value only when the condition (A12) holds. Un- 
der the condition (A12), we only have to consider three 
cases (i) |A^| = |A^,|-i ^ 1, (ii) r = r', |A^| = 1, and (in) 
T ^ r', |A^| = |A^,| = 1. When jA^ | > 1 and |A^J < 1, 



7- 



IT 



because finite values of 7^ and 7^^ cause 
the divergence of \d^\ at j = 00 and j = —00, respec- 
tively. Thus, 7^*7^ in Eq. (A6) must be zero in case 
(i). In case (iii), (A^) and (A^,) are accidentally degen- 
erate. Since this degeneracy is lifted by an infinitesimal 
change in the energy E, case (iii) can be excluded. Since 
?7^/ ^ + rj'^ ^, is nonzero only in case (ii), Eq. (A5) is rep- 
resented by 

ATM 

lUj) = J2^MZ';.r)b!^\'+lMZ'^r.,-r)b'-r\' ■ (A14) 



Equations (A5) and (A6) do not depend on j because 
A^,*A^ = A^*A^ = |A^|2 = 1 in case (ii). 

When the vectors fj^ are normaUzed as Ini(Z[T-| i^-i) — 1 
and Ini(Z_|7-| _|t-|) = —1, we can show 



N' 



N^ 



hot = Y.(Wr? - \lLr?) = E(l7?l' - \^^r?) 



N'- 



Y.{\i??-\i%?) 



(A15) 




inner tube outer lube 

I I 




attached probe; 

FIG. 1. TDWNT 



Here, conservation of probabihty guarantees that Itot is 
common to regions /, Z?, and O. Eq. (A15) is equivalent 
to *5^^*5^^ = 1, *5^o*5^0 ^ 1 and *^^o*^/o ^ i_ 



APPENDIX B: RELATION BETWEEN THE 
NORMALIZATION AND GROUP VELOCITY 



From Eqs. (13) and (21), we can obtain 
Differentiating Eq. (Bl) with respect to k, 



(Bl) 



^ dk dk ^ dk 

•^" dk ^^^ ^" dk 



(B2) 



Using Eqs. (12), (13), (21) and (42), Eq. (B2) is repre- 
sented by 



" dk ^^ 



2Re 



dff^-^ „ dB^'-^' 

E!f 



dk 



dk 



g'r 



Aj = l ^ 



(B3) 



by 



From Eqs. (51), (A7) and (A8), Eq. (B3) is represented 



dE^^ 



(B4) 




FIG. 2. Structure [fi] superimposed on its initial structure. 
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With Eqs. (23), (36) and (B4), Eq. (30) can be repre- 
sented by 



FIG. 3. The dispersion lines near the neutral Fermi level 
Ep are numbered based on the index r under the conditions 

Ki ((-l)^lf ) = ^^M (-(-l)^l^) - 4°' and r'-§. > 0. 



%0{k?rpOJ{k',)gJ{ki) 



\lm{ZO^)lm{Zl^)\a^Ak^ 



(B5) 



When the normalization |lm {Z^^) | = |lm (Z^ ^) | = 1 
is used. 



%0^k?rpO.'{ki)gJiki) 



a^ko - kiy 



(B6) 
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FIG. 4. The dispersion lines without interlayer interaction 
for M = 1 (open circles) and Al = 2 (closed diamonds), where 
M denotes the integer parameter of Eq. (12). 
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FIG. 5. Effective Hamiltonian elements defined by Eq. (14). 
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FIG. 7. The dispersion lines for (a) Ae = -0.015 (b) Ae = 
and (c) Ae = 0.020 Ry. At fc = 0.657r/a, e_i = -0.022 + Ae 
and e2 ~ —0.028 + Ae. The crosses represent k^{E) and 
kl-{E) without interlayer interaction. The dispersion lines 
with interlayer interaction are shown by circles for the ex- 
act k^ {E) and by solid lines for the approximate k^' (E). 
The wavenumber k is calculated for the discrete energies 
E = -0.428 + 0.004J where j = 0, 1, 2, ■ • • , 14. 
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FIG. 6. Schema for the matrix P. 
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tlie number of unil cells in the overlap region (i) 
FIG. 8. The transmission rate Ti,i at an energy E = —0.4 
Ry. The values of Ae are the same as in Fig. 7. The 
circles represent the exact T\^\ with Eqs. (46) and (47). 
The crosses represent the exact Ti,i when go a-nd go in Eqs. 
(46) and (47) are replaced by Q^ ■ The average defined by 
r''"'^(L) = (r(L - 1) + T{L) + T{L + l))/3 are shown by 
the dot-dashed and dashed lines for the circles and crosses, 
respectively. The approximate formula (38) is shown by thick 
solid lines. The thin solid represent Eq.(38) where EB '^ in 
Eq.(23) is replaced by zero, i.e., the interlayer overlap matrix 
is set to zero. The thin solid lines are shown only for the first 
period. 




4 10 20 30 40 50 
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FIG. 9. The same calculation as in Fig. 8 for the transmis- 
sion rate T2,2- The crosses and dashed lines are omitted. The 
triangles show Eq. (56). 
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